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We propose a measure of order in the context of nonequilibrium field theory and argue that this 
measure, which we call relative configurational entropy (RCE), may be used to quantify the emer- 
gence of coherent low-entropy configurations, such as time-dependent or time-independent topologi- 
cal and nontopological spatially-extended structures. As an illustration, we investigate the nonequi- 
librium dynamics of spontaneous symmetry-breaking in three spatial dimensions. In particular, 
we focus on a model where a real scalar field, prepared initially in a symmetric thermal state, 
is quenched to a broken-symmetric state. For a certain range of initial temperatures, spatially- 
localized, long-lived structures known as oscillons emerge in synchrony and remain until the field 
reaches equilibrium again. We show that the RCE correlates with the number-density of oscillons, 
thus offering a quantitative measure of the emergence of nonperturbative spatiotemporal patterns 
that can be generalized to a variety of physical systems. 



I. INTRODUCTION 



Spatially-localized solutions to partial differential equations pQ play a central role in many current models of 
particle physics [5] and condensed-matter physics [5] . Some of the most studied of these solutions are topological and 
nontopological defects that may appear during spontaneous symmetry breaking 4J. As is well-known, topological 
defects owe their stability to the nontrivial topology of the vacuum manifold [5], while nontopological defects owe 
theirs to the conservation of a global charge [6]. In both cases, solutions involve one or more interacting fields and 
are either time- independent or have a simple harmonic time-dependence ~ expfiwt], as is the case of Q-balls [7] and 
other nontopological solutions [8]. 

Another class of localized-energy solutions, known as oscillons, exhibits nontrivial time dependence During 
the past decade or so , oscillons were shown to exist in models with a single self- interacting scalar field in dimensions 
d < 6 [T2] and in many models with gauge fields, Abelian [13] and nonAbelian [T41 [15], including the Standard 
Model [16] . They can be exceedingly long-lived, self-supported by their nontrivial interactions due to feedback from 
parametric resonance |15L I17j . Oscillons can also play an important role in cosmology, as reported in recent work 
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Both topological and nontopological defects and oscillons can be thought of as attractors in field configuration 
space |21j : given that certain dynamical constraints are satisfied, for a broad range of initial conditions the system 
will evolve toward these solutions. This is not necessarily surprising, since these configurations are solutions to the 
equations of motion. Semi-classically, we expect them to dominate the path integral. 

We may, however, ask whether there is another way to quantify the existence of nontrivial spatially-localized 
solutions in field configuration space. Can we think of spatially-coherent field configurations as ordered states in an 
informational sense, in analogy with the Shannon entropy of information theory 22, 23]? That is, given the set of 
field modes that are allowed by the constraints of the model (i.e. initial and boundary conditions), do they carry a 
special informational signature that can be quantified? 

To this end, we recently proposed a measure of configurational entropy, based on the Fourier transform <j){k) of 
square- integrable, bound functions <f>(x) (24] ■ Leaving the details for Section III, here it suffices to say that with this 
configurational entropy, we can establish a correlation between the energy of a localized-field configuration and its 
associated configurational entropy. In particular, we can show that departures from the solution of the eom will have 
correspondingly larger relative configurational entropies. 

In the present work, we take this approach one step further, applying it to nonequilibrium fields. Our recent 
treatment of Ref. [54] was for static solutions to the field equations, such as Id kinks and 3d bounces. Here, we will 
add time dependence, computing the relative configurational entropy during spontaneous symmetry breaking in the 
context of a 3d scalar field model. "Relative" refers to a comparison between the entropy of the field at some time and 
the entropy of the initial state, which we take to be a thermal state at temperature T. We will compute the change 
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in relative configurational entropy as the field is tossed out of equilibrium during the symmetry-breaking process. 
Previous work has shown how, for certain types of quench, oscillons naturally emerge during symmetry breaking 
[T5"l ITT] . Our results show quantitatively that the emergence of these localized coherent structures coincide with the 
largest departures from equilibrium and that they carry the most information content, in the sense defined in Section 
III. We are thus proposing a measure of order in field configuration space. Furthermore, we are able to establish a 
direct correlation between the emergence of ordered spatiotemporal structures, in effect the number density of such 
structures, and the relative configurational entropy. 

The work is organized as follows. In Section II we describe the model and its lattice implementation. In Section III we 
define the relative configurational entropy (RCE) and apply it to spontaneous symmetry breaking. We describe how the 
RCE provides a measure of ordering in the system as compared to the maximally-disordered initial thermal state and 
how this order correlates with the existence of spatially-coherent structures. In Section IV we present our conclusions 
and plans for future work. In the Appendix, we provide some technical details of the lattice implementation. 



II. THE MODEL 



We consider a (3+l)-dimensional scalar field theory with Lagrangian density 
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and the tree level potential given by 



C=\{d^f-V{ct>), (1) 



V{4>) = ^<t> 2 - f ^ + (2) 

where the parameters m, a, A are positive-definite and temperature independent. We use h = c = kg = 1 and rescaled 
variables <f)' = (j)v\/m, x' M = x^m, a' = a/(mvA) to write the potential as V(4>) = (m 4 / \)V(<j)'), with 



W) = 




We henceforth drop the primes and work with the rescaled variables. In this work we only consider the values a = 
and a = 3/2. The first case corresponds to a potential with a single minimum at </> = and the second describes 
a double-well potential with degenerate minima at <fi — and (j) — 2. For a — the potential is .Z2 -symmetric. 
For any other value of a this symmetry is broken. In the context of 2d [17] and 3d [15] models, it has been shown 
that when quenching a thermalized field from the symmetric to the double-well potential (here, from a = to 
a — 3/2), large-amplitude fluctuations about the vacuum state give rise to oscillon formation as the system evolves 
towards thermal equilibrium in the new potential. In other words, coherent, spatially-extended configurations develop 
spontaneously when an initially featureless system is tossed out of equilibrium. The mechanism behind this process is 
well understood: coherent oscillations of the field's zero mode parametrically amplify higher fc-modes and the resulting 
energy transfer triggers the formation of oscillons [111 H7J ■ This remains true when the expansion of the universe is 
incorporated into the dynamics jT5] EH] • We note that we could have investigated the traditional symmetry-breaking 
mechanism with a double- well symmetric about = 0. In this case, the Z2 symmetry-breaking would lead to spinodal 
decomposition and the formation of domain walls. 



A. Lattice Implementation 

We simulate the formation of oscillons using a cubic lattice with iV 3 = 256 3 points, periodic boundary conditions, 
lattice spacing dx = 0.5 and time step dt = 0.01. We prepare the initial thermal state using standard Langevin 
dynamics [55] 

+ 70-V 2 = -|^+C, (4) 

where V((j)) is the rescaled potential in Eq. [3] with a — and £ is a Markovian noise with two-point correlation 
function obeying the fluctuation-dissipation relation 
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(C(x, i)C(x', f )) = 2 7 T<5(x - x')<5(i - 0, (5) 

where T is the temperature parameter characterizing the initial state and we take 7 = 1. The system was evolved 
until thermalization was achieved, indicated by the onset of equipartition with every mode having average kinetic 
energy T/2, after which the potential was quenched from a = to a = 3/2 and the field's coupling to the heat bath 
was removed (7 — > 0). So, after the quench the dynamics is conservative. The evolution of the field was done using 
a symplectic velocity Verlet algorithm and the Laplacian was discretized with a second-order accurate, forth-order 
isotropic stencil using all 26 neighbors of a 3 x 3 x 3 cube around a point [26 . We have checked that our results 
showed no particular dependence on the choice of box size, lattice spacing and timestep. 



B. Hartree Approximation and Oscillon Emergence 



Thermal fluctuations in 4> will change the potential V(<p) which, to leading order in perturbation theory, can be 
approximated by the Homogeneous Hartree Approximation |27j . Since the Hartree approximation assumes that the 
fluctuations of the field remain Gaussian throughout its evolution, it works well just before and after the quench for 
all temperatures. For low temperatures, it remains valid at all times. We can thus derive the Hartree potential by 
writing the field as <j> — <fi + S<fi and then averaging over all fluctuations 6<j> to get Vjj — (V(4> + 6<f)). Under the 
Hartree assumptions we have (6<f>) — and (5(f) 2 ) = /3, where (S(f> 2 ) is the mean square variance of the field and is 
proportional to the temperature parameter T. Suppressing the bar, the Hartree potential becomes 



V H (<l>) = -aP<t>+(l + ~l3^-a^- + ^.. (6) 

When thermal equilibrium is reached through the coupling to the heat bath, the field modes in momentum space 
satisfy 



(l'Mfc)l 2 ) = ^2~T~33T2~' ( 7 ) 



•h 

with the Hartree mass given by m 2 H = V^(0) = (l + |/3). 

On the lattice, Eq. [7] has to be adjusted for lattice effects due to the cf> dependence on the discretization and the 
inherent UV cutoff on the lattice. For the lattice UV cutoff ir/dx, f3 can be analytically obtained in terms of the 
temperature T as [15] 
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^^dx < 8 > 

The continuous dispersion relation u> 2 = k 2 + m 2 H has to be modified to take into account the field dependence on 
the discretization scheme. For the isotropic discretization we use here, the radially-averaged dispersion relation will 
be given by uj 2 — k 2 s + m 2 Hl with [5S] 



c 4 2 r /2 r /2 

k 2 a = — — -k \ cos9 / \2c3 [cos [k cos 4> cos 9dx] + coslk sine/) cos 9dx] + coslk sin 9dx]] 

dx z Trdx z Jo Jo 

+4c2 [cos [k cos (f> cos 9dx] cos [k sin (j> cos 9dx] + cos [k cos <fi cos 9dx] cos [k sin 9dx] + cos [k sin cos 9dx] cos [k sin 9dx] ] 
+8ci cos [k cos cj> cos 9dx] cos[fc sin <j> cos 9dx] cos[k sin 9dx]\ d<j)d9, (9) 

and c% = 1/30, C2 = 1/10, C3 = 7/15, C4 = —64/15 being the discretization coefficients of the Laplacian. With these 
lattice effects taken into account, the two-point correlation function at equilibrium becomes 



(l<O fc )i 2 > = 



T 

l + 9T/(8Trdx)' 



(10) 
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FIG. 1: Radially-averaged two-point correlation functions for the field at T 
data at equilibrium and the analytical approximation (^^"(fc)] 2 } (solid line) of Eq. 

along with the two-point correlation function data at time 44m -1 after the quench (diamonds). Modes with k < 1.5m get 
significantly amplified while the rest remain in thermal equilibrium throughout the simulation. 



0.25. On the left, we show the simulation 
On the right, the same data is plotted 



In the left part of Fig. [T]we show the radially-averaged two-point correlation function for the field (j> after it reaches 
thermal equilibrium at temperature T = 0.25. Squares denote the data from the numerical simulation, and the black 
solid line is the theoretically predicted spectrum (|</4q tt (fc)| 2 ), adjusted for lattice effects. The averaging for low fc-modes 
is less accurate because, on the lattice, there are fewer modes to compute the average of the power spectrum. As k gets 
larger, there are more modes with the same value of k and the agreement with the theoretically-predicted spectrum 
is evident. On the right we plot the equilibrium spectrum (squares) and the spectrum at t = 44m -1 (diamonds) 
after the quench. Modes with k < 1.5m have been parametrically amplified, whereas modes with k > 1.5m have 
remained in thermal equilibrium. The low fc, out of equilibrium modes that are parametrically amplified after the 
quench are the ones responsible for oscillon formation. Oscillons emerge synchronously after the quench as the system 
is tossed out of equilibrium, and then slowly disappear as the field evolves towards its new equilibrium state [T3J [T7] . 
A visualization of the process can be seen at [25] . In Fig. [2} we show a snapshot of the field at t = AAm 1 after a 
thermal Wiener filter has been applied to remove the high fc-modes still in equilibrium, with oscillons appearing as 
localized spikes in the field. 



III. INFORMATION CONTENT OF COHERENT FIELD CONFIGURATIONS 

Reverting back to Fig. [I] we observe a marked difference in the two-point correlation function of field modes between 
the equilibrium (left) and nonequilibrium phases (right). In particular, as we emphasized before, low fc-modes with 
k < 1.5m are greatly amplified when the field is away from equilibrium |15j . In this section, we propose a measure 
to quantify this nonequilibrium amplification of low fc-modes. We further show that this measure, having a natural 
interpretation as an entropy in field configuration space, correlates with the number of coherent field configurations 
(oscillons, in the case studied here) that emerge as the system is tossed out of equilibrium. Thus, we propose that our 
entropic measure can be used to quantify the emergence of complexity in field theory, if by complexity we understand 
the appearance of spatially-localized coherent field configurations, in contrast with the structureless (disorganized) 
thermal state. 

Following our recent work |24) . we define the modal fraction of a field 0(x, t) in Fourier space at time t as 
Equivalently, in equilibrium we have 



(12) 
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FIG. 2: Snapshot of the field cj> at t = A4m~ x after the system is quenched from a single to a symmetric double-well. The 
simulation size shown here is L 3 = 128 3 . The field was initially thermalized at T = 0.25. Oscillons (spatial ordering) appear as 
spikes about the zero mode of the field. Three isosurfaces are shown at <j> = 0.5, 1.3, 1.6 in purple, cyan and red, respectively. 
As the system evolves towards its new equilibrium, the spatiotemporal ordering is gradually lost and oscillons subsequently 
disappear. A full visualization of the process can be seen at [29], where it is also clear that for early times oscillons emerge in 
synchrony (time ordering). 



with |</() eq (k)| 2 given by Eq. [7] Note that with this definition, g(k) depends on the temperature T only through the 
Hartree mass correction term /?. In the lattice implementation, g(k) is computed with ^^(k)] 2 as defined in Eq. 
|10| The linear dependence on the temperature is cancelled by our choice of normalization, which gives the modal 
fraction units of m~ 3 . Next, we define the dimcnsionless relative configurational entropy (RCE) of the power spectra 
as 



S f (t) = J S f (k,t)d 3 k, (13) 
where the relative configurational entropy density «Sf(k, t) is given by 



5 / (k,t) = /(k,<)ln^^. (14) 

Our definition of the RCE is a field-theory version of the Kullback-Leibler divergence (KLd) commonly used in infor- 
mation theory to compare two probability distributions P and Q of a discrete random variable: Z?kl = Y)Pj ^°-(Pi/Qi) 
[30] . In information theory, the KLd gives a measure of the expected number of extra bits required to code samples 
from P using a code based on Q. Usually, P represents the "true" data or a precisely computed distribution, while 
Q represents a theory, model, or approximation of P. 

By comparing the modal fraction of the field's Fourier transform at time t with that of the thermal state, we 
obtain a measure of the amplification of the low fc-modes responsible for oscillon formation, a "distance" in Fourier 
configuration space from the thermal state. (In principle, other mode expansions could be used, although, as we 
argued in Ref. [24j . for localized fields the Fourier transform is the most natural.) Modes that remain in equilibrium 
throughout the evolution of the field (here, roughly for \k\ > 1.5m) have f{k,t) = g(k) and Sf(k,t) = 0, and thus 
do not contribute to the RCE. The larger the modal fraction f(k,t), the larger its associated relative entropy density 
S f (k,t). 

The RCE thus provides a clear measure of the departure from equilibrium. Furthermore, it peaks where the field 
is most organized into coherent spatial structures. This is consistent with the notion that the farther a system is 
from equilibrium, the farther it is from satisfying equipartition. In effect, since the thermal state has maximum 
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entropy and hence no information (all modes have the same average energy-as equipartition determines) the RCE 
gives a measure of information in field configuration space: peaks in the RCE correspond to peaks in information-rich 
coherent structures. In the case here of a scalar field with an attractive self-interaction, the attractor point of its 
dynamics, the oscillon, is the farthest that it can be from equipartition. This explains why the RCE peaks when 
oscillons are present, with the amplitude of the peaks correlating directly with the number of oscillons, as we will see 
below. 




FIG. 3: Relative entropy density S/(|fe|, i) and relative entropy Sf(t) as a function of time for a simulation of initial temperature 
T = 0.25. Initially the relative entropy is zero, corresponding to the time just after the quench when the system is still in 
equilibrium. As parametric resonance takes place, spikes in the low k part of the relative entropy density emerge, signaling 
the formation of oscillons. The amplitude of the spikes correlate with the number of oscillons formed, as explained below. As 
the system evolves back to equilibrium, the oscillons disappear and the relative configurational entropy goes back to zero. The 
wave vector magnitude k has units of m and time has units of m _1 . 



In Figj3]we plot the RCE density (left) and RCE (right) as a function of time for the same simulation with T = 0.25 
of Figs. Ill and ^) Since the field starts in thermal equilibrium, the RCE density, Sf(\k\,t), is zero everywhere initially. 
After the quench, energy is transferred to low-fc modes and spikes in Sf(\k\,t) begin to appear. These are clearly seen 
after integrating over the k modes, as shown in the plot of the RCE on the right. 

In order to substantiate the claim that the RCE gives a quantitative measure of the emergence of spatiotemporal 
structure in the system, in Fig. El we plot both Sf(i) and the number density of oscillons n osc (t) — N osc (t)/V. It 
is quite clear that the spikes in Sf(t) coincide with the synchronous emergence of oscillons, and that the higher the 
amplitude of Sf(t) the larger the number density of oscillons present. (An advantage of defining the relative entropy 
to be dimensionless is that it can be consistently used for similar simulations in lattices of different size. When the 
simulation size is increased while working with the same initial temperature, the power spectrum will look the same 
but the number of oscillons formed will be proportionally larger. Therefore, our definition of the relative entropy is a 
consistent measure of the number density of oscillons formed.) 

As an illustration, the highest peak in the two-point correlation function, shown in the right-side plot of Fig. [T] 
appears at time t — AAmT 1 . The isosurface snapshot of Fig. [2] shows the richness of spatial structure in the field at 
that time. This is also the time when the highest peak in the RCE appears, and corresponds to the maximum in the 
number density of oscillons, as is clear from Fig. [4] Eventually, as the system evolves towards its final equilibrium 
state, the relative entropy goes slowly back to zero, signaling the disappearance of coherent structures in the field. 



For low temperatures, the Hartree approximation is valid for all times and Eq. [10] describes well the equilibrium 

as 



spectrum. For larger temperatures (T > 0.30), the Hartree approximation breaks down and we cannot use Eq. 10 



descriptive of the equilibrium state of the system. Although it will still be a good base for studying the emergence of 
structure in the field, it should be used with care. 

To verify that the correspondence between the RCE and the number density of oscillons holds for the temperature 
range where oscillons appear in the system, we extract the maxima of the RCE and the corresponding oscillon number 
density at that time for a range of initial temperatures 0.10 < T < 0.29. (For example, for T = 0.25 this would be 
at t — 'AAm~ 1 in the simulation displayed in Fig. [4]) For each value of T we perform an ensemble average over 15 
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FIG. 4: Relative configurational entropy Sf(t) (continuous line) and number density of oscillons n OS c(t) (dashed line) as a 
function of time for a simulation of initial temperature T — 0.25. Initially, both the relative entropy and the number density 
of oscillons are zero, since the system starts in equilibrium. After the quench, a clear correlation is seen between the spikes in 
Sf(t) and the maxima of n OBC . The wave vector magnitude |fe| has units of in and time has units of m _1 . 



simulations and plot the results in Fig. [5j The vertical axis on the left displays the ensemble-averaged maximum 
value of the RCE, S'max, with the data shown in circles, while the right vertical axis displays the corresponding 
ensemble-averaged oscillon number density n osc with the data depicted by diamonds. The error bars show the 
standard deviation of the ensemble. For low temperatures 0.10 < T < 0.15, the fluctuations on the field do not have 
large enough amplitude to lead to oscillon formation and the RCE is zero. As the temperature increases, parametric 
amplification of the oscillon-related k modes triggers the formation of oscillons. The relationship between the RCE 
and the number density of oscillons is evident: higher temperature leads to both higher RCE and larger oscillon 
number-density. For temperatures T ~ 0.28 — 0.29 we begin to see a discrepancy between the two values, the reason 
being that the Hartree approximation starts to diverge from the true equilibrium spectrum. (For reference, symmetry 
restoration occurs roughly at T ~ 0.32.) We note that the results of Fig. [5] are independent of the simulation volume. 



IV. CONCLUSION AND DISCUSSIONS 



We investigated the noncquilibrium dynamics of symmetry breaking in the context of a 3d real scalar field model with 
a double-well potential. Preparing the system in a parity-invariant initial thermal state, we break the Z2 symmetry 
by adding a cubic term to the potential. As a result, coherent spatiotemporal structures emerge, an ensemble of long- 
lived oscillons. We proposed a measure to quantify the emergence of spatiotemporal order, which essentially counts 
the modes out of equilibrium. This measure, which we called relative configurational entropy following our previous 
work of Ref . |24j , provides an accurate description both of the departure from equilibrium and of the emergence of 
coherent structures in the field. "Relative" here refers to computing an entropic distance in field configuration space 
from a baseline which we took to be the thermal state. This way, we are able to provide the informational content of 
nonequilibrium field structures, in particular of coherent states that emerge during spontaneous symmetry breaking. 
We could just as easily have studied the information content of domain-wall formation had we used a different initial 
state centered at the maximum of the double-well potential. 

The measure we proposed here should be generalizable to many different contexts. For example, it should be 
possible to apply it to models with gauge fields, or in nonrelativistic applications of interest to condensed matter 
physics, such as Ginzburg-Landau models of superfluids and superconductors. Furthermore, a similar measure should 
also be of interest in cosmological contexts where topological and nontopological structures appear due to symmetry 
breaking driven by the cosmic expansion. Finally, it would be interesting to see if such ideas could be extended 
to compute the relative configurational entropy of metric spaces, thus providing a possible measure of gravitational 
entropy. Work along these lines is currently in progress. 
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FIG. 5: Ensemble-averaged values of the maxima of the RCE and of the oscillon number density for a range of different initial 
temperatures. For low temperatures T < 0.15, no oscillon formation is possible and the RCE is zero. As the temperature 
increases, the maxima of the RCE increase in amplitude at the same rate as the maxima of the oscillon number density. For 
temperatures T ~ 0.28 — 0.29 the rates begin to diverge as the Hartree approximation is no longer valid. 
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V. APPENDIX: ANALYTICAL ESTIMATE OF THE MODAL FRACTION NORMALIZATION AT 

EQUILIBRIUM 

The normalization of the modal fraction at equilibrium in Eq. [IT] can be computed analytically. This will provide 
an estimate of how important lattice effects are for different lattice parameters. We assume spherical symmetry in 
momentum space and impose a UV cutoff A to perform the integral 

/ \(t>^(k)\ 2 d 3 k= [ T 2 4nk 2 dk = 4ttT (a - m H tan' 1 —\ , (15) 
J Jo k + m H \ ran ) 

where T is the temperature and ran is the Hartree mass. The value of the UV cutoff A is chosen appropriately 
to describe the discretization employed in the numerical calculations. Here, we pick the value A = y/in/dx which 
corresponds to the largest wavevector represented in a cubic lattice of lattice spacing dx. Geometrically, this choice 
makes the integration volume of Eq. [15] to be the smallest sphere that fully covers the cubic lattice in momentum 
space. 

We evaluate Eq. [l5]for T = 0.25 and several choices of lattice spacing dx. The results are shown in Table [IJ For 
large values of dx, the analytical approximation is in good agreement with the numerically calculation; but as the 
lattice spacing is decreased, the two values begin to diverge. To see why the analytical approximation underestimates 
the value of the integral, note that k 2 grows faster than the numerically-calculated k 2 s of Eq. M and this becomes 
increasingly more noticeable for smaller lattice spacings and hence larger volume of integration. This further justifies 
our incorporation of lattice correction effects in the dispersion relation in Eq. |10[ as the discrepancy is of the order of 
10% for dx — 0.5 used throughout this paper. 
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dx 


f\^\k)\ 2 d 3 k 


f \<j>™(k)\ 2 d 3 k 


0.9 


14.8874 


14.3853 


0.7 


20.9985 


19.6331 


0.5 


32.0204 


29.1701 


0.3 


57.611 


51.5837 


0.1 


184.715 


164.262 



TABLE I: Comparison between the numerical and analytical values of the modal fraction normalization at equilibrium for 
different values of lattice spacing. The analytical approximation given by Eq. [15] works well as long as the lattice spacing is 
not very small. For smaller values of dx, the lattice effects on the dispersion relation calculated in Eq. [9] have to be taken into 
account. The lattice spacing dx is given in units of mT x and the normalization factors have units of [Tin]. 
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